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Abstract We show how to detect entanglement with 
criteria built from simple two-body correlation terms. 
Since many natural Hamiltonians are sums of such cor- 
relation terms, our ideas can be used to detect entangle- 
ment by energy measurement. Our criteria can straight- 
forwardly be applied for detecting different forms of mul- 
tipartite entanglement in familiar spin models in thermal 
equilibrium. 



1 Introduction 

Entanglement is an important non-classical phenomenon 
in quantum mechanics which plays also a crucial role in 
the novel field of Quantum Information Theory. While 
for pure quantum states it is equivalent to correlations, 
for mixed states the two notions differ. In this general 
case, an iV-qubit quantum state is entangled if its density 
matrix cannot be written as a convex sum of product 
states 



i 



(i) 



States of the form Eq. are called separable. Based on 
this definition, several sufficient conditions for entangle- 
ment have been developed. In special cases, e.g. for 2x2 
(two-qubit) and 2x3 (qubit-qutrit) bipartite systems 
E] and for bipartite multi-mode Gaussian states [3] even 
necessary and sufficient conditions are known. 

However, in an experimental situation usually only 
limited information about the quantum state is avail- 
able. In this case, only those approaches for entangle- 
ment detection can be applied which require the mea- 
surement of not too many observables. One of such ap- 
proaches is using entanglement witnesses. They are en- 
tanglement conditions which are linear in expectation 
values of observables. The theory of entanglement wit- 
nesses has recently been rapidly developing [2E]|S1|H]. It 
has been shown how to generate entanglement witnesses 



that detect states close to a given one, even if it is mixed 
or a bound entangled state [7||H] • It is also known how to 
optimize a witness operator in order to detect the most 
entangled states [§]. Apart from constructing witnesses, 
it is also important to find a way to measure them. Opti- 
mal measurement of witnesses have been studied in Refs. 
|1UU11II12| . Recently, witnesses have been developed to 
detect entanglement in physical systems in the thermo- 
dynamical limit |13II14II15| . 

Entanglement witnesses can not only be used to de- 
tect entanglement experimentally, but can also be used 
to characterize the entanglement of a multipartite quan- 
tum state. As we will see later, in the multipartite set- 
ting several different classes of entanglement occur, and 
entanglement witnesses can be used to decide in which 
class a given state is [7]. 

In this paper we ask what we can do for systems of 
very many particles, e.g., for spin chains in thermal equi- 
librium. Entanglement in spin chains has already been 
extensively studied Hi . In Section 2 we discuss how to 
detect entanglement in general in spin models based on 
the ideas presented in Ref. 13 . In Section 3 we study 
the detection of different types of multipartite entan- 
glement in these systems as discussed in Ref. F° r 
this aim, we determine what the important questions are 
from this point of view in spin systems in the thermo- 
dynamical limit. Then we look for appropriate entangle- 
ment witnesses, which are easy to construct and study 
multipartite entanglement with them 17 . 

2 Bipartite entanglement 

Let us consider first the two-qubit case. The simplest 
expression which can be used for entanglement detection 
must contain at least two correlation terms 



W:=AWAW+BMBP 



(2) 



where A^ and Bk are operators acting on qubits /c = 1,2. 
For simplicity, let us consider A^ and Bk with eigenval- 
ues ±1. Now, if we want to use W for entanglement 
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detection, we have to make sure that 

M(W) 9 < inf (W) , (3) 

where the right hand side of the equation is minimized 
over the set of product states V. (Minimization over all 
mixed separable states would lead to the same value due 
to the convexity of the set of separable states.) Eq. J2J 
expresses the fact that the minimum of (W) must be 
larger for separable states than for quantum states in 
general. For that it is necessary to have 

[A k ,B k ]^0, (4) 

for k — 1,2. Here [..] denotes the commutator. Eq. J3J 
expresses the fact that we have to measure two different 
observables at each party. For entanglement detection in 
an experiment, the ratio of the two minima in Eq. J3J 
must be the largest possible. It is straightforward to see 
that this is the case if we choose operators such that 

{A k ,B k } = 0, (5) 

for k — 1,2. Here {..} denotes the anticommutator. An 
example for such an operator is then 

hxY:=xWxW+YWY<?\ (6) 

where X and Y denote Pauli spin matrices. For this 
operator the minimum of the expectation value is 

M(h XY h = -2- (7) 

The state giving the minimum is the two-qubit singlet 

h/, s ):=i=(|10>-|01». ( 8 ) 

For this state (X^X^) = (yWyC 2 )) = -1. The mini- 
mum for product states can be obtained as follows. For 
product states we have 

(h XY ) = (XW)(X( 2 )) + (Y«)(Y( 2 )) > -1, (9) 

where the last inequality follows from the Cauchy-Schwarz 
inequality and knowing that {X^) 2 +{Y^) 2 < 1. Among 
operators with three correlation terms used for entangle- 
ment detection the following form is optimal 

h H := X«Jf( 2 ' + y«y( 2 ) + z^z {2 \ (10) 

The minimum of the expectation value of hn is —3. For 
separable states the minimum is — 1 which can be proved 
similarly as it has been done for Hxy- 

Now let us move to the A-qubit case. Consider the 
expression 

N 

H XY := J^lWl(W)+y ( % (w) , (11) 

k=l 

where J > is the coupling constant and according to 
the usual assumption for a periodic boundary condition 



qubit (N+ 1) is identical to qubit (1). This is the Hamil- 
tonian for the XY chain. The minimum for separable 
states is now 

inf (Hxy}* = -JN. (12) 

This comes from knowing that for product states each 
term in the summation in Eq. (jTT|) is bounded by — 1 
as we have seen it before. The minimum for quantum 
states can be obtained from numerical calculations since 
the XY model is solvable Similarly, we can define 

N 

H H ■= J^2x {k) X {k+1) + Y {k) Y {k+1) + Z {k) Z {k+1 \ 
fe=l 

(13) 

This is the Hamiltonian for the Heisenberg chain. The 
minimum for separable states is 

inf (H H } = - JN. (14) 

The minimum for quantum states can be obtained for 
large N as [T3] 

mf(H H )v =-4^1n2- ^ J NJ w —1.77NJ. (15) 

For our spin chain Hamiltonians the ratio between 
the minimum for general quantum states and the mini- 
mum for separable states is smaller than for Eqs. (|6I10|) 
since there is not a quantum state saturating all two- 
body correlation terms. In fact, it is easy to see that 
there is not a Hamiltonian built from two-body corre- 
lations such that its unique ground state saturates all 
correlation terms and this ground state is true multipar- 
tite entangled |2"U] . 

What are the advantages of the expressions Eq. (|11|) 
and Eqs. ijT3|) in detecting entanglement? They are easily 
measurable locally, since they are the sum of only a few 
two-body correlation terms. Moreover, in some physical 
systems Eq. (I13f) can directly be measured as the aver- 
age nearest-neighbor correlation, or as the energy of the 
system if this system can be described by a Heisenberg 
Hamiltonian. 

The previous ideas can straightforwardly be applied 
to spin chains in thermal equilibrium |13) . Let us con- 
sider the Heisenberg Hamiltonian in an external mag- 
netic field H HB ■= H H + BJ2 k Z (k) . For this Hamil- 
tonian, it is easy to bound the minimum for separable 
states ^31- Any time (Hh) is below this value, we know 
that the thermal state is entangled. This is demonstrated 
in Fig.^ It shows the nearest-neighbor entanglement vs. 
B and the temperature T. The entanglement of forma- 
tion was computed from the concurrence EH- Light color 
indicates the region where the thermal ground state is 
detected as entangled based on the ideas discussed be- 
fore. As one can see, there are regions with Ep > 
which are not detected. However, it can be seen in the 
Figure that when the system contains at least a small 
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Fig. 1 Heisenberg chain of 10 spins. Nearest-neighbor entan- 
glement as a function of magnetic field B and temperature 
T. 

amount of entanglement (~ 0.07) the state is detected as 
entangled. Note that the sharp decrease of the nearest- 
neighbor entanglement around B cr i t = 4 for T = is 
due to a quantum phase transition. 

3 Multipartite entanglement 

In this section we will discuss how to detect multi-party 
entanglement by measuring the operators described be- 
fore. Our motivation is that entanglement for many par- 
ticles is qualitatively different from the two-party case, 
and many new phenomena arise |22|. There are several 
possibilities to classify entanglement of many parties. We 
are now looking for the terminology which is appropriate 
for spin chains of many particles. 

Let us first recall the notion of genuine multipartite 
entanglement. A pure state \ip) of a quantum system of 
N parties is called fully separable if it is a product state 
for all parties, 1^) = \4>i) <8 |<fe) <8> ■■■ <8> \4>n) ■ It is called 
biseparable, when a partition of the iV parties into two 
groups A and B can be found, such that the state is a 
product state with respect to this partition, namely 

IV) = \4>a) ® l<M • (16) 

If this is not the case, the state is called genuine mul- 
tipartite entangled. Note that the vectors \4>a) and \4>b) 
are allowed to contain entanglement within their parti- 
tion. Thus, to prove genuine multipartite entanglement, 
it does not suffice to exclude full separability. 

For mixed states, these definitions can, as usual, be 
extended via convex combinations. Indeed, the defini- 
tion of full separability was already given in Eq. 1(1)1. 
A mixed state is biseparable, whenever we can write 
g = J2iPi\' l Pi)(' l l ; i\ with biseparable \tpi) and some prob- 
abilities pi. Here, the states are allowed to be bisep- 
arable with respect to different partitions. 

Another approach to classify multipartite entangle- 
ment asks whether multipartite entanglement is neces- 
sary to form a given state JS] . In this approach, a state 



(a) . . . 

(b) . . . 




Fig. 2 (a) Chain of two-qubit singlets, (b) The same state 
shifted by one qubit to the right. The mixture of these two 
states is two-producible, that is, does not need three-qubit 
entanglement when produced from pure states by mixing. 

producible by k-party entanglement (or fc-producible, 
in short) if we can write the state as a tensor product 

\ip) = <g> \<f> 2 ) ® ... ® \<t> m ) , (17) 

where the states \<pi) are states on maximally fc-qubits. 
In this definition, a two-producible state does not con- 
tain any multipartite entanglement, since it suffices to 
generate the two-qubit states \4>i) to arrive at the state 
. In addition, we say that a state contains genuine 
k-party entanglement if it is not producible by (k — 1)- 
party entanglement. This definition can be extended to 
mixed states as before via convex combinations. Again, a 
mixed state which is fc-producible requires only the gen- 
eration of fc-party pure entangled states and mixing for 
its production (see also Fig. 2). Consequently, a mixed 
state q contains fc-party entanglement, iff the correla- 
tions cannot be explained by assuming the presence of 
(k — l)-party entanglement only in the pure subsensem- 
bles. 

The notions of genuine multipartite entanglement and 
producibility are not completely independent. For exam- 
ple, the states containing A-party entanglement are just 
the genuine multipartite entangled states and the one- 
producible states are fully separable. If one can show 
that a reduced state of k + 1 qubits is genuine multi- 
partite entangled, then this implies that the total state 
is not fc-producible, while the converse is in general not 
true. 

For spin chains of macroscopic size, it is in general 
very difficult to prove that the total state is genuine 
A-partite entangled via energy measurements. This is 
due to the fact that the notion of genuine A-partite en- 
tanglement is extremely sensitive to the properties of a 
single qubit. Indeed, in order to prove genuine multi- 
partite entanglement, one has to exclude the possibility, 
that one single qubit can be separated from the remain- 
ing A — 1 qubits. However, multipartite entanglement in 
the reduced states of small numbers of qubits can eas- 
ily be detected, as we will see. Moreover, if the reduced 
state is multipartite entangled then the state is not two- 
producible. 

Now let us see our results for the XY model and the 
Heisenberg chain. The proofs for the following theorems 
are given in Ref . |15| . We always assume periodic bound- 
ary conditions and that the number of spins A is even. 

Theorem 1. Let g be an A qubit state whose dy- 
namics is governed by the XY-Hamiltonian in Eq. 1111) . 
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If g is one-producible (fully separable), then 

(H XY ) > ~JN (18) 

holds. If (Hxy) < —JN this implies that there are two 
neighboring qubits such that their reduced state is en- 
tangled. For two-producible states 

(Hxy) > -pN (19) 

holds. If (Hxy) < —9/8JN the state contains thus tri- 
partite entanglement and if 

1 H - /2~ 

(Hxy) < ^— JN ~ - 1-207 JJV (20) 

then there exist three neighboring qubits i, i+ l,i + 2 

such that the reduced state gi,-<+i,i+2 of these qubits is 

genuine tripartite entangled. 

For the Heisenberg model, we can state the following: 
Theorem 2. Let g be an N qubit state with the 

Heisenberg Hamiltonian of Eq. If g is one-producible 

(fully separable), then 

(H H ) > -JN (21) 
holds, while for two-producible states 

(H H ) >-\jN (22) 

holds. Thus, if (Hh) < 3N/2 the state contains genuine 
tripartite entanglement. Furthermore, if 

(H H ) < - JN w -1M8JN (23) 

then there are three neighboring qubits such that their 
reduced state is genuine tripartite entangled. 

Let us see an example. Consider the state shown in 
Figure Eta) 

|*») = \A) ® IV> S > <8 \ip s ) ® (24) 

where the two-qubit singlet is defined in Eq. JBJl. It is 
easy to see that state |<P S ) saturates the inequality Eq. (|2"2"|l 
It is not surprising, since it is a two-producible state. Let 
us now define operator S which shifts the qubits by one, 
i.e., 

S |ai) <g> |a 2 ) <8> |a 2 ) <8> ... <8> |ajv) 

= |ajv) <g> |ai) <g> |a 2 ) ® ... <8> |ajv-i) ■ (25) 

Consider the state 

p m :=~(\$s)(<i> s \+S\<I> s )($ s \S^. (26) 

This state is the mixture of the singlet chains depicted 
in Figure EI a ) an d (b). p m is not fully separable and 
is not a product of single-qubit and two-qubit density 
matrices. Moreover, the state p m has a negative par- 
tial transpose with respect to each partition. However, 



Table 1 Threshold temperatures Tc2,T«3 and Tc3 for a 
Heisenberg chain. The parameters as set to J = fcs = 1. See 
text for details. 



N 


2 


4 


6 


8 


10 


Tc'2 


7.28 


3.45 


3.21 


3.18 


3.18 


Tc3 




2.10 


1.75 


1.65 


1.62 


Tb3 




1.85 


1.46 


1.32 


1.26 



p m also saturates the inequality Eq. 1)2 2H and it is also 
two-producible. That is, three-qubit entanglement is not 
needed to create it, and it contains no multipartite en- 
tanglement. 

The previous results can straightforwardly be used 
for obtaining a limit temperature for the different forms 
of entanglement. We define thus the temperatures Tr 2 , 
Tr3i Tc2 and Tc3 below which either reduced states of 
two or three parties are entangled or the total state con- 
tains two or three-party entanglement. Obviously, T# 2 = 
Tc2 > Tc3 > Trs has to hold here. These temperature 
bounds are shown for a Heisenberg chains of a couple of 
spins in Table 1. As expected, the values for Tq% = Tin 
coincide with the ones of Ref. ^2] ■ The given values for 
Tc3 and Tbz show that in the Heisenberg chain of ten 
spins at fc^T sa J multipartite entanglement plays a role, 
namely at least one reduced state is genuine tripartite 
entangled and the total state contains tripartite entan- 
glement. 

4 Conclusions 

We discussed how to construct entanglement conditions 
using two-body correlations. This implies, that typical 
Hamiltonians as appearing in the XY model or the Heisen- 
berg model can serve for entanglement detection in spin 
models. Also different forms of multipartite entangle- 
ment can be detected in this way. 

A natural continuation of our work lies in the exten- 
sion of our bounds to other systems. Here, spin systems 
in two or three dimensions as well as frustrated systems 
are of interest. Furthermore, it would be also desirable 
to derive energy bounds also for higher classes of multi- 
partite entanglement, e.g. three-producible states. 
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